import numpy as np
import numpy.matlib as ml

n = 80
pmax = 20
w = np.sqrt(list(range(pmax)))

n = 5

b = 25.867325

A = np.matrix([
  [ 1.075334, 0.526197, -3.608734, 0.657207, 0.990262 ],
  [ 0.526197, -0.867184, 3.377548, 3.454253, 1.561950 ],
  [ -3.608734, 3.377548, 1.450808, 1.426643, 1.431982 ],
  [ 0.657207, 3.454253, 1.426643, 2.818069, 3.047428 ],
  [ 0.990262, 1.561950, 1.431982, 3.047428, 0.977788 ],
]);

C = ml.eye(n);
